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Abstract. The transition between the Main Sequence and the Red Giant Branch in low mass 
stars is powered by the onset of the CNO burning, whose bottleneck is the ''^NCp,/)''©. The 
LUNA collaboration has recently improved the low energy measurements of the cross sec- 
tion of this key reaction. We analyse the impact of the revised reaction rate on the estimate 
of the Globular Clusters ages, as derived from the tumoff luminosity. We found that the age 
of the oldest Globulars should be increased by about 0.7-1 Gyr with respect to the current 
estimates. 
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1. Introduction 

Globular Clusters (GCs) represent the oldest resolved stellar populations. Their age practically 
coincides with the time elapsed since the epoch of the formation of the first stars in the Universe 
and provides an independent check of the reliability of standard (and non standard) cosmological 
models. Moreover, the age spread in the GC system is a primary indicator of the time scale of the 
halo formation. Among the various methods to date stellar Clusters, the most reliable and widely 
adopted is that based on the measurement of the luminosity of the turnoff (i.e. the bluest point 
on the main sequence). This dating technique requires the knowledge of the Cluster distance, 
the light extinction along the line of sight and the chemical composition (see |Gratton et al, 2003| 
for an exhaustive analysis of the present state of the art). In addition, a reliable theoretical cal- 
ibration of the turnoff luminosity-age relation (TOL-A) is needed. This relies on our knowl- 
edge of the physical processes of energy generation (e.g. nuclear reactions) and transport (e.g. 
opacity) taking place in H burning low mass stars. An adequate description of the thermody- 
namics of stellar matter is also required. Finally we have to consider any mechanism capable 
to modify the internal chemical stratification (once again nuclear reactions, convective mixing, 
rotational induced mixing, microscopic diffusion or levitation induced by radiation pressure). 
[Chaboyer et al (1996)| discuss the influence of various theoretical uncertainties on the calibration 
of the turnoff' luminosity-age relation and conclude that the total uncertainty due to the theory 
may be confined within 0.5 Gyr. 

This paper is devoted to the evaluation of the impact on the theoretical calibration of the 
Globular Clusters ages of the improved determination of the rate of the key reaction ^'^N(p,y)'^0, 
as obtained by the LUNA collaboration l |Formicola et al, 2003bt . Since thermonuclear reactions 
are responsible for chemical modifications occurring in stellar interiors and supply most of the 
energy irradiated from the stellar surface, the estimated stellar lifetime depends on accurate mea- 
surements of their rates. In the last few years, many efforts have been spent in improving these 
measurements at energy as close as possible to the Gamow peak, namely the relevant energies at 
which nuclear reactions take place in stars. This is a mandatory requirement for the calibration 
of stellar ages. 

The main sequence stars presently observed in Globular Clusters have masses smaller than 
that of the Sun. As in the Sun, these low mass stars burn H in the center, mainly through the pp 
chain. However, towards the end of their life, when the central hydrogen mass fraction becomes 
smaller than about 0.1, the nuclear energy released by the H burning becomes insufficient and 
the stellar core must contract to extract some energy from its own gravitational field. Then, the 
central temperature (and the density) increases and the H burning switches from the pp chain to 
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the more efficient CNO cycle. Thus, the departure from the Main Sequence is powered by the 
CNO cycle, whose bottleneck is the '"^N(p,y)'^0 reaction. The luminosity of the turnoff depends 
on the rate of this key reaction: the larger the rate, the fainter the turnoff. On the contrary, the 
total lifetime is only marginally affected by a change in the CNO, because it is mainly determined 
by the rate of the 'H(p,e^, Vf)'H. As a consequence, an increase of the CNO rate would imply 
fainter turnoff points, for a given age, or younger ages, for a given turnoff luminosity (see also 
[Chaboyer et al, 1998| l. Note that an equivalent effect may be caused by the enhancement of the 
CNO abundances ( |Rood, 198l||Salaris, Chieffi & Straniero, 1993> . 

In the next section we remind the new measurements of the stellar cross section of the 
''*N(p,y)'^0 reaction. Then, in section 3 we present the revised turnoff luminosity-age relation 
(TOL-A). We show that this revision leads to systematically larger estimates of the age of the 
Globular Clusters. Implications for cosmology are briefly discussed in the conclusive section. 

2. The updated ^"^NCp.y)!^© reaction rate 

The minimum energy explored in nuclear physics laboratories before LUNA was ~240 keV, 
which is well above the range of interest for the stellar CNO burning (~20-80 keV). Therefore, 
the reaction rate used in stellar model computations is largely extrapolated, in a region where the 
resonant structure of the '^O compound nucleus is particularly complex. The rates reported by 
the popular compilations ( |Caughlan & Fowler, 19881 CF88, and [Angulo et al, 1999| NACRE), 
which are based on the cross section measurements obtained by [Schroder et al (1987)| are very 
similar. In particular, the astrophysical factor at zero energy is S(0) = 3.2 + 0.8 keV b (NACRE). 
The main contributions to S(0) come from the transitions to the groundstate in '^O and to the 
subthreshold state at Ecm = -504 keV. It is the existence of this subthreshold resonance that 
makes the extrapolation very uncertain. Recently [Angulo & Descouvement (200 1)| re-analyzing 
the Schroder's experimental data by means of a R-matrix model, report a significant lower S(0), 
namely 1.77 + 0.20 keV b. The main discrepancy concerns the contribution of the captures to 
the '^O groundstate, which has been found 19 times smaller than the value quoted by Schroder 
and adopted by NACRE and CF88. We emphasize that the large discrepancy among different 
analyses based on the same data set is a clear demonstration of the inadequacy of the low energy 
extrapolation for this reaction. 

The LUNA collaboration' has significantly improved the low energy measurements of this re- 
action rate ( |Formicola et al, 2003b| l. We used a 400 keV facility ( |Formicola et al, 2003at , which 
is particularly well suited when reaction y-ray lines up to ^ 7.5 MeV have to be measured 
with very low intensities. Cosmic background is strongly suppressed by the mountain shield- 
ing and low intrinsic activity detectors are employed. The explored energy window ranges 
from 390 keV down to 135 keV, i.e., significantly closer to the astrophysical relevant energy 

' LUNA is an acronym of the Laboratory for Underground Nuclear Astrophysics, operating at the LNGS 
of Assergi, Italy. 
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Fig. 1. Comparison between evolutionary sequences obtained with different rates of ''*N(p,'y)'^0: 
NACRE (dotted line), CF88 (solid line) and LUNA (dashed Une). The stellar mass is 0.8 Mq and 
the metalhcity is Z=0.0003 

than any previous experiment. The fit of the new data by means of a R-matrix model leads to 
S(0) - 1.7 + 0.1 (stat) + 0.2 (sys) keV b. In the following we use this result to revise the 
calibration of the turnoff luminosity-age relation. 

3. Globular Clusters ages 

New stellar models have been computed with the same code described in 
[Straniero, Chieffi & Limongi (1997)| but updating the rate of ''*N(p,y)'^0. We recall that 
this code includes an improved equation of state (which properly takes into account the degree 
of degeneracy of the electrons and the electrostatic interactions), the most recent compilation 
of opacity for stellar interiors l |Iglesias, Rogers & Wilson, 1992| [Alexander & Ferguson, 1994| l 
and microscopic diffusion ( |Thoul, Bahcall & Loub, 19941 1. Fi gure 1 shows an example of the 
evolutionary track obtained by adopting different rates for the ''*N(p,')')'^0 reaction. The tracks 
obtained with the CF88 rate practically coincide with the one obtained with the NACRE rate. 
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Fig. 2. Isochrones for Globular Clusters obtained with different rates of the ^"*N(p,7)'^0 reaction: 
CF88 (solid) and LUNA (dashed). The brightest isochrone (of each set) is the youngest (10 Gyr), 
while the fainter is the oldest (18 Gyr). The age step between two adjacent isochrones is 1 Gyr. 
The metallicity is Z=0.0003 (or [M/H]=-1.82). 

On the contrary, the tumoff and the subgiant branch of the sequences obtained by adopting the 
new rate are substantially brighter. 

Isochrones have been computed for two sets of stellar models, the first based on the old 
CF88 rate and the second based on the revised LUNA rate. We have explored the whole range 
of chemical composition covered by the galactic GC system. In particular, the mass fraction of 
metals (the metallicity) has been varied between Z = 0.0001 and Z = 0.006, which corresponds 
to [M/H]= -2.3 and [M/H]= -0.5 ^. Some examples of the comparison between old and new 
isochrones are shown in Figure 2 and 3. As expected, the lower rate of ^^N(p,y)^^0 leads to 
brighter and bluer tumofF points (for a given age). When a given tumofF luminosity is considered, 
the revised isochrones imply systematically older ages, namely between 0.7 and 1 Gyr. 

To compare our isochrones to the available photometric studies of globular cluster stars, we 
have transformed luminosities and efl'ective temperatures into magnitudes and colors by means 

^ standard spectroscopic notation. 
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Fig. 3. Isochrones for Globular Clusters obtained with different rates of the ''*N(p,y) '''O reaction: 
CF88 (sohd) and LUNA (dashed). The brightest isochrone (of each set) is the youngest (10 Gyr), 
while the fainter is the oldest (18 Gyr). The age step between two adjacent isochrones is 1 Gyr. 
The metallicity is Z=0.003 ([M/H]=-0.52) in the right panel. 

of model atmospheres provided by Castelli et al. (1997). The accuracy of the new isochrones 
in reproducing the morphology of the observed color-magnitude diagrams, has been checked by 
selecting two clusters which are representative of the oldest component of the galactic halo. 
The first test is illustrated in Figure 4. Isochrones for Z - 0.0003 ([M/H]= -1.8) and age 
13, 14 and 15 Gyr are over imposed to the color magnitude diagram of NGC 6397. A similar 
test, but for NGC 5904 (M5), is reported in Figure 5, where the isochrones have Z = 0.001 
([M/H]= -1.3). Photometric data are from the ground-based database published by Rosenberg 
et al. (2000). In both cases, the new isochrones match the overall color-magnitude diagram 
well at 14 Gyr, with a bona fide uncertainty of +1 Gyr. Similar results were obtained by 
[Straniero, Chiefli & Limongi (1997)1 with the old (CF88) isochrones, but, in that case, the best 
reproduction of the observed diagrams required 13 Gyr (see their Figure 1 1). 

The following relation for the Globular Cluster age (fg in Gyr), as a function of the V magni- 
tude of the turnoff point (M™) and the metallicity ([M/H]), has been derived. 
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Fig. 4. Test to the CMD of the metal-poor cluster NGC 6397. The new isochrones with 13, 14 
and 15 Gyr are reported. Their metalUcity is Z=0.0003 ([M/H]=-1.8). We adopt (m-M)v=12.58 
and E(B-V)=0. 18. The data are from Rosenberg et al. (2000). 



logfg = -1.0146- 0.273 1[M/H] + 0.03032[M/H]M(,'' 

The standard deviation of the estimated age is A log tg 
ages ranging between 10 and 18 Gyr and [M/H] between 

4. Implications for cosmology 

The recent developments of accurate measurements of the fundamental cosmological parameters 
allow us to derive a very precise age of the Universe: to = 13.7 + 0.2 Gyr ( |Spergel et al, 2003) 1. 
This result has been obtained in the framework of a ACDM model and it is based on the mea- 
sures of three fundamental parameters: Hq, whose best determination has been obtained by the 
Key HST Project (|Freedman et al, 2001| l, Q, measured by WMAP ( jSpergel et al, 2003> , and the 
ratio ^m/^a, constrained by the observation of type la supernovae in high redshift Galaxies 
( jPerlmutter et al, 1 999| [Smith et al, 1998} . The galaxy clustering shape measurements also con- 



- 0.00058([M/H]M™)2 + 0.4801*?^;^^ 

= 0.005. This relation can be used for 
-2.3 and -0.5. 
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Fig. 5. Test to the CMD of the intermediate metallicity cluster NGC 5904. The new isochrones 
with 13, 14 and 15 Gyr are reported. Their metallicity is Z=0.001 ([M/H]=-l). We adopt (m- 
M)i/=14.41 and E(B-V)=0.02. The data are from Rosenberg et al. (2000). 

strain Om ( |Percival et al, 200"!) . It is obvious that any systematic uncertainty affecting just one 
of these experiments would imply a revision of this estimate of the age of the Universe. For 
example, it has been argued that the light curve of an SNe la might depend on the chemical 
composition and/or the mass of the progenitor star. In this case, the commonly assumed simi- 
larity between nearby and high redshift SNe la, in spite of the different stellar populations of 
their host Galaxies, could have induced a systematic error in the evaluation of D.m/^a (see e.g. 
IDormnguez, Hoflich & Straniero, 2001) 

In this context, an independent determination of the age of the Universe, may (or may not) 
confirm the standard cosmological model that emerges from the experimental cosmology. At 
present, the most reliable dating technique is the one based on the TOL-A relation for the oldest 
stellar systems of the Milky Way, the Globular Clusters. The standard cosmological model also 
predicts that the H reionization, which should coincide with the epoch of the first star formation, 
occurred between 100 and 400 Myr after the Big Bang (95% CL, |Spergel etal, 2003> . Such a 
delay must be also considered. 
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An exhaustive comparison between stellar and cosmological ages requires a detailed statis- 
tical analysis taking into account all sources of errors (experimental and theoretical). This is 
beyond the purpose of the present paper and will be presented elsewhere. Let us limit our discus- 
sion to the expected implication of the revised ''*N(p,y)'^0 reaction rate. We have shown that the 
revised ages of the Globular Clusters are older, about 0.7-1 Gyr, than those previously claimed. 
It is worth to note that, in the framework of the ACDM model, an equivalent increase of to might 
be obtained by reducing Ho (~5%) or Q.m (~8%). These variations are, in any case, within the 
experimental errors. 

Gratton et al (2003), by means of the TOL-A relation derived from models taking into ac- 
count the effect of microscopic diffusion, but computed adopting the old '"^N(p,y)'^0, conclude 
that the age of the oldest Galactic Clusters is 13.4 Gyr (+0.8 random, +0.6 systematic). When 
the age increment implied by the revision of the ''*N(p,y)'^0 is barely added, the best fit to the 
age of the oldest GCs should increase above 14 Gyr. This revised lower limit of the age of the 
Universe strengthens the need of a positive cosmological constant. In the case of a flat Universe 
(D. = 1) and assuming Ho = 0.72 Km s ' Mpc~' jFreedman et al, 2001t , a Universe older than 
14 Gyr would imply Q.m < 0.22 or, adopting an uncertainty of 1 .4 Gyr in fo ( |Gratton et al, 2003) , 
D.M < 0.35. Note that this upper limit for the matter density is independent of the SNe la obser- 
vations. Alternatively, by coupling our result with that of the high redshift SNe la, we may relax 
the assumption on the geometry of the Universe to derive a stringent constraint for the Hubble 
constant. Indeed, taking Hoto = 0.96 + 0.04 ( |Tonry et al, 2003^ , the present lower limit for to 
would imply Hp < 67 Km s ' Mpc ' (or Ho < 74 within 1 cr in fo), in good agreement with 
72 + 8 obtained by the Key HST Project ( |Freedman et al, 200 !> . 
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